Nonequilibrium Josephson-like effects in wide mesoscopic S— N— S junctions 
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Mesoscopic superconducting-normal-metal-superconducting (S-N-S) junctions with a large separation 
between the superconducting electrodes (i.e. wide junctions) exhibit nonequilibrium supercurrents, even 
at temperatures for which the equilibrium Josephson effect is exponentially small. The second harmonic 
of the Josephson frequency dominates these currents, as observed in recent experiments. A simple 
description of these effects, in the spirit of the Resistively-Shunted-Junction model, is suggested here. 
It is used to calculate dc I—V characteristics, and to examine the effects of various types of noise and 
of external microwave radiation (Shapiro steps). It is found that the nonequilibrium supercurrents are 
excited when the junction is driven by a dc bias or an ac bias, or even by external noise. In the case of 
junctions which are also long in the direction perpendicular to the current flow, thermodynamic phase 
fluctuations (thermal noise) alone can drive the quasiparticles out of local equilibrium. Magnetic flux 
is then predicted to be trapped in units of $o/2 = hc/Ae. 

PACS numbers: 74.50.+r, 74.40. +k, 74.80.Fp, 73.23.Ps 



I. INTRODUCTION 

Wide mesoscopic superconducting-normal-metal- 
superconducting (S-N-S) junctions exhibit a spectrum 
of low-lying electronic "Andreev bound states" |}J , which 
depends strongly on the phase-difference <j> between the 
two superconductors ||. This </>-dependence persists 
when the temperature T is raised, even if the normal- 
metal coherence length, £n, becomes much smaller than 
the distance L x between the two superconducting elec- 
trodes. In this high-temperature regime the equilibrium 
Josephson coupling is exponentially weak, but the effects 
of quantum-mechanical coherence of the electrons decay 
only as 1/T (Ref. As the spectrum depends on <j>, 
and hence on time t (the derivative dfi/dt is proportional 
to the voltage V), non-equilibrium (NEQ) effects occur, 
including supercurrents which are the topic of this work. 

NEQ proximity effects have been studied intensively 
for more than 20 years [0, but experiments in the 
mesoscopic regime have only recently become possible. 
The latter include Refs. which displayed intriguing 
frequency-doubling effects and thus motivated the devel- 
opment of the theory presented here (an interpretation 
in terms of the frequency-doubling effects predicted in 
the early eighties by Spivak and coworkers j^] was ini- 
tially proposed, but had to be rejected Theoreti- 
cally, NEQ phenomena have been studied in both clean 
0-[t3| and dirty p|- |i"s[] l imits. The approach taken 
here, following e.g. Ref. is to consider the junction 
as a mesoscopic system, or "electron box". The spe- 
cial superconducting "walls" of this box impose a time- 
dependent boundary condition — the external phase- 
difference 4>. Once the ^-dependence of the spectrum of 
the system is known, the supercurrents are simply given 
by the "slopes" of the occupied levels, dE n /d<p. In fact, 
it is shown that, for evaluating the NEQ supercurrent, 
all of the different levels may be replaced by a single rep- 




FIG. ffl: Energies of Andreev bound states E n , equally 
spaced in n, for a diffusive S-N-S junction (inset, Ec — 
TiD/L x <C A), plotted as a function of the superconduct- 
ing phase-difference <p. A phase-dependent mini-gap is 
proximity-induced at the Fermi level (e = 0). Blowup 
(schematic): the avoided crossings of individual disorder 
realizations, which are not resolved here as the spectrum 
is dense. 

resentative state, at a relatively low cost in terms of loss 
of accuracy (~15%). This results in a particularly trans- 
parent description of short junctions, L y ~ L x (see inset 
in Fig. |l|), which may be used in a wide variety of appli- 
cations, in analogy to the well-known resistively shunted 
junction (RSJ) model. The simplicity of the model al- 
lows us to consider for the first time NEQ-dominated 
long junctions, L y 3> L x . Although a different language 
is used, the agreement with the quasiclassical Green's 
functions approach used in Refs. jff],[l4|,[l6) is complete, 
as we will show explicitly. 

We will assume that the normal-metal part of the 
junction is mesoscopic in the sense that the dephasing 
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length lu, is much larger than L x , but that the dephas- 
ing time t v is sufficiently short to prevent the distinction 
between individual quantum levels, H/t v 5, where S is 
the single-particle level spacing (in the present context, 
a system with resolved levels, h/r^ ~ 5, would best be 
called microscopic) . This implies that we are discussing a 
two- or three-dimensional system, in which the number 
of transverse channels N± ~ Uvf/L x 5 is very large (vf 
is the Fermi velocity). The effects of electron-electron 
and electron-phonon interactions in the normal-metal 
part will be ignored pofl , except for their inclusion in 
the dephasing rate 1/t v and the energy-relaxation rate 
1/te, which are treated as phenomelogical parameters. 
However, the BCS interaction is assumed to dominate in 
the superconducting electrodes, so that they form good 
reservoirs for Cooper pairs, with well-defined phases (or 
at least phase difference, (f>) and a large gap A for single 
particle excitations. 

We focus on the "high-temperature regime", Eq <C 
r« A, where Eq is the Thouless energy or correlation 
length of the spectrum (the temperature T is in energy 
units, /cb = 1). The first inequality here is equivalent to 
£n C-ti, and implies that the equilibrium Josephson ef- 
fect is exponentially supressed. The second implies that 
single-particle excitations are "bound" to the normal- 
metal region and do not occur in the superconducting 
electrodes — any electron or hole which approaches the 
N-S interface is either reflected normally or Andreev re- 
flected back into the normal-metal region (Andreev re- 
flections are coherent processes by which an electron can 
add a Cooper pair to the superconductor while evolving 
into a hole, or vice-versa, a hole can remove a Cooper 
pair and evolve into an electron; these processes con- 
vert current carried by Cooper pairs to current carried by 
single-particle excitations; re-diagonalizing the Hamilto- 
nian in their presence, i.e. solving the Bogoliubov — de 
Gennes equations, gives the Andreev bound states, which 
form the single-particle excitation spectrum of the sys- 
tem pi]). The Thouless energy Eq is defined as Ti over 
the time it takes an electron to cross the junction. In the 
dirty limit, l c \ <C L x with l e \ the elastic mean-free-path, 
the Thouless energy is given by Eq = hD/L x , with D 
the diffusion constant. In this case, the normal-metal 
coherence length is £n = \JTiD/2hT. We will also dis- 
cuss the clean or quasiballistic case, Z c i 3> L x , where 
Eq = hv-p/Lx and £n = fivp jlixT . However, we will 
not discuss the super-clean limit, Z c i 3> N±L X , for which 
fr/Vei ~ 5 where t c \ — I c \/vf — in that case disorder is 
irrelevant and the spectrum may become seperablc. Our 
discussion will also be limited to slow energy-relaxation 
rates, 7i/te <C Eq. This inequality is already implied 
by our previous assumptions, because te > t v and the 
dephasing length is large, l v ^> L x . 

The layout of this paper is as follows. The model is de- 
veloped in Sec. II, including simple applications to NEQ 



ac and dc Josephson-like effects. Further applications, 
including long junctions, are analyzed in Sec. III. The 
model is adjusted to describe a clean, quasiballistic sys- 
tem in Sec. IV, and conclusions are briefly discussed in 
Sec. V. The consistency of the present treatment with the 
standard but more opaque Green's function approach is 
demonstrated in the appendix. 

II. EQUILIBRIUM AND NONEQUILIBRIUM 
SUPERCURRENTS 

A mesoscopic junction may be characterized by a den- 
sity of states, v{t,<j)), which depends on the excitation 
energy e (measured from the Fermi surface) and the su- 
perconducting phase-difference </>. In this section we will 
take this spectrum as input, and calculate the resulting 
supercurrents. 

Let us begin by displaying the spectrum for a simple 
diffusive junction, for specificity. The disorder-averaged 
density of states, v(e , </>), can be calculated from the Us- 
adel equations It is convenient to define the en- 

ergy E n {<f)) of the nth level: n — f^ n v(e,(f)) de, with 
£L n = — E n , see Fig. |l| (negative excitation energies are 
included here for clarity of presentation; they represent 
the "hole branch" of the excitation spectrum). The cal- 
culation assumes time-reversal symmetry, with the pair 
potential vanishing in N (no electron interactions), and 
perfect N-S interfaces (i.e. no potential barriers or Fermi 
velocity mismatch). A mini-gap is induced in N by the 
proximity of S. Its size is E g ~ 3.1Eq at <f> = 0, and it 
closes at <fr — tt and reopens periodically Jl| . We are con- 
sidering only the disorder-averaged spectrum, because 
the fluctuations in E g and E n are only ~ <5, and are asso- 
ciated with currents that are smaller than the ones con- 
sidered here by a factor of 5/Eq. This, the ^-periodicity, 
and the scale Eq of the structure in the spectrum near 
the Fermi level, are generic to all junctions to be con- 
sidered here. Other features, such as the existence of a 
strict gap, are not generic and will not be relied upon — 
for example, due to the great sensitivity of the spectrum 
to various types of symmetry breaking, adding spin-flip 
scattering with h/r s f ~ Eq can close the mini-gap. An- 
other concrete example — a clean junction — will be 
considered in Sec. IV. 

The current / in the junction is related to the spec- 
trum. In the presence of a voltage V = (Ti/2e)(d(f>/dt) 
(the Josephson relation), the energy of the electronic 
system, J2 n E n fn, is time-dependent Dividing the 
power by the voltage, one finds / = /s+-^n with 

2e ^ dE n ^ E n (df n \ /1N 

n=— oo n— — oo 

where Ig is independent of V, and will be called the su- 
percurrent, and In is proportional to V (dissipative), 
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and will be called the normal current. As is by now 
well-known (see e.g. Ref. the same distinction be- 

tween a dissipative current and a persistent current can 
be made in a normal-metal mesoscopic ring with a non- 
trivial boundary condition <f> imposed by an Aharonov- 
Bohm flux. In that case the contributions of the differ- 
ent levels E n ((j)) to the persistent current Is do not add 
up in a systematic manner, the resulting currents are 
much smaller than in S-N-S junctions, and in a disorder- 
averaged model of the type considered here, would vanish 
altogether. 

As we are considering a system with a dense spectrum, 
the occupations /„ do not evolve adiabatically (in con- 
trast to Ref. (ill), but diffuse in energy 0] with a co- 
efficient D E (e,<t)) ~ G N V 2 S, where G N ~ (e 2 /h)E c /S 
is the conductance of the normal-metal part in the ab- 
sence of proximity effects. This can be argued quali- 
tatively by observing that in a time t, a quasiparticle 
crosses the junction ~ tEc/h times, and undergoes An- 
dreev reflections which shift its energy by ±eV every time 
it encounters the N-S interface, implying difusion with 
L>e ~ { e V) 2 D / L 2 X - For T > E c regions of weak prox- 
imity effects, |e| > Ec, dominate in the normal current, 
and we will thus set In = G^V. 

While the occupations /„ follow the slope of the E n (</>) 
curves and diffuse in energy, they are also subject to 
electron-electron and electron-phonon interactions. We 
use the relaxation time approximation p4| , 

^ - --(/„-«) , (2) 
at te 

where / cq (e) = 1/ (l + exp(e/T)) , and consider only 
small voltages, where the term D E (d 2 f/de 2 ) which was 
dropped from the right hand side of Eq. (^|) can indeed 
be neglected, and heating is unimportant (T is time- 
independent). This implies te-De <C E c , or eV <C 
y/ Ech/rE, which is more stringent than the condition 
eV <C Ec required for the spectral density itself to be 
a meaningful quantity (at higher voltages the spatial de- 
pendencies of the relevant Green's functions must also be 
taken into account). 

In equilibrium, I s = de j(e, 0)/(l+exp(e/T)) = 
I cq (4>), where j(e) is the "Josephson current density", 
j(E ni cf>) on v(E n ,(p){dE n /d<j)), which is analytic in the 
upper half of the complex e plane (a positive imagi- 
nary part of e corresponds to a dephasing rate, h/r v , 
which would smooth out any singularities). By contour 
integration one finds the well-known Matsubara sum: 
I oq = 2mTY^ =l j{iu m ), where u> m = (2m-l)nT. 
At high temperatures, even the smallest Matsubara fre- 
quency is ^> Ec and has a decay time fr/wi < L 2 /D, 
yielding an exponentially small j(iuj m ). For physical 
quantities of this form, thermal averaging is thus "equiv- 
alent" to dephasing. 

In NEQ situations contour integration cannot be used, 
and the currents are not exponentially small. The large- 
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FIG. p|: Supercurrents for a dc voltage bias: (a) Is vs. t, 
and (b) amplitudes Ik of its harmonic decomposition at 
t 3> te (with k > 0; phases not shown), for five dc volt- 
ages labeled by v = {2e/K)VT E . The Josephson period is 
ij = h/2eV (for I ncq see below). 

T behavior of / cq implies oscillations of j(e) such that 
/ de e k j(e, <fi) = for any odd power k, as can be seen 
by expanding / eq in powers of 1/T. In the diffusive case, 
j(iuj) decays exponentially with ^Juj/Ec on the imag- 
inary axis, and correspondingly j(e) oscillates and de- 
cays rapidly with ^/~eJ~Ec (the curves of Fig. |l| change 
their character repeatedly at higher e, occasionally hav- 
ing shallow maxima at cf> — tt rather than minima). As 
can be seen from Eq. the deviations of /„ from equi- 
librium change sign in rhyme with these oscillations (and 
decay with e), and thus the integrand of the NEQ part 
of the supercurrent, 2 J Q dej (f n — feq), does not change 
sign and cannot be affected by any cancellations. 

When a dc voltage is applied to the junction, dej)/dt = 
2eV/H = const., Eqs. (|l|) and (||) give rise to ac super- 
currents, see Fig. ||. Here T Ec was used to take 
/cq(e) = I — e /4T, and J eq = 0. The second harmonic 
dominates because Is involves products of dE n /d<f> and 
(fn—feq), an d both factors oscillate at the Josephson fre- 
quency. This "frequency doubling" is an essential feature 
in identifying these effects experimentally ||. 

The limits of small and large voltages relative to 
the inelastic scatterring rate may be analyzed analyt- 
ically p]Jl6l|. For small (2e/K)VTE, the deviations of 
the occupations from equilibrium are linear in the volt- 
age, /„ ~ / eq - TE,[2e/h)V(df e Jd4>), where df cq /dcf> ~ 
(-l/AT){dE n /d4>) for T > E n . The NEQ supercurrent 
is then disipative — the energy is dissipated by the fast 
ralaxation te and not stored in the system. From Eq. (Q), 
one finds that 
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with / = G(0)y. Despite the fact that Is is dissipative in 
this limit, we still use Eq. ([!]) to distinguish the compo- 
nent of the current we call "normal" from the "supercur- 
rent" or "persistent current" which remains finite even 
when the voltage vanishes and equilibrium is restored, 
provided the temperature is sufficiently small. 

In the opposite limit of large dc voltages, {2e/h)Vr-E ^> 
1, the individual occupation probabilities /„ relax to the 
time-average of their equilibrium value, 



fn 



1 

2^ 



E n 

4T ' 



(4) 



where E n is the phase-averaged value of E n (<f). The 
supercurrent is then obtained directly from Eq. (|l|), and 
is "purely reactive" . In fact, the rate of dissipation in 
this limit is proportional to 1/te and may be calculated 
from Eq. (J2J) . As it is independent of the voltage, the dc 
component of the current dacays as 1/V at large voltages 
(within the domain of validity of the model) see Fig. ||. 

The presence of the sum over n (the energy integration) 
in Eq. (Q) makes calculations according to this model rel- 
atively cumbersome. In fact, it may often be superfluous, 
because the relaxation-time approximation used here al- 
ready ignores possible e-dependencics, e.g. in te. One 
may also note that only the shape of the E n {4>) curves 
is significant — a simple shift in energy does not affect 
their contribution. We will thus adopt a much simpler 
scheme, in which only the first N states above (or be- 
low) the Fermi level are considered to have a significant 
^-dependence, and are all represented by a single curve, 
E Iep (<f>): 

I S ~ 7 oq (0) + 2^N^ (/ - f cq (E Icp )) , (5) 

where f = jt Yln=i fn ^ s ^ ne avera g e occupation of the 
represented band of N levels |25|. For the system of 
Fig. 0, E lcp = 3.4E C \/l + 0.7cos(V) and N = 10E C /S 
(including a factor of 2 for spin) reproduce the result of 
the full sum over n very well — deviations are less than 
15%, for 85% of the harmonic amplitudes shown in Fig. || 
(the results plotted are from the simpler model). 

Often, I(t) rather than V(t) is known. In considering 
the conventional Josephson effect in such situations, it is 
customary to use the resistively shunted junction (RSJ) 
model pq] , which describes the dynamics of </> through 
(h/2e)(d(j)/dt) = (I— Is)/Gn, with Is = I c sin</> as appro- 
priate for the equilibrium case. This corresponds to over- 
damped motion of a "representative point" in a "tilted 
washboard potential" , F(4>) = — (h/2e)(I c cos (f)+I4>). To 
generalize this to NEQ situations where the occupation 
/ needs to be followed, one replaces the expression for I s 
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FIG. |: The I—V curves and differential resistances (up- 
per insets) for a dc current bias in the NEQ model, for 
four values of the relaxation time te; compared to the 
RSJ model with I c = I ncq , and to I = G^V (dotted 
lines). Lower inset: The effective free-energy function 
F(<f),f), for / = 0, is a "winding valley" with minimum 
at / = E((f>) and a parabolic cross-section. 

by Eq. (@), and describes the dynamics of / by Eq. (Q). 
This is associated with the free-energy landscape 



F(<f>,f) - E neci (E - ff - (h/2e)I<j> 



(6) 



(see inset in Fig. y), where the overdamped motion is 
given by 



h d(f> 
2e~dl 



1 2edF 
Gn h dej) 



df_ 

dt 



1 



dF 



2TE-Eneq df 



(J) 



Here E and / are rescaled quantities, E lcp ((f>) = AE + B 
and / = | — (Af + B)/4T, with the requirement -1 < 
E < 1 fixing A and B (in the present case, A ~ l.SEc). 
The current scale is given by / ncq = (2e/h)NA 2 /4T, and 
the energy scale is -E neq = (h/2e)I ncq ~ A.IE^/ST, only 
a factor of Eq/3.5T smaller than the T = Josephson 
coupling energy, Ej(0) ~ 1AEq/S. 

The dc I — V curves of this model, obtained by taking 
I(t) = const., are displayed in Fig. |[ using units of 7 ncq , 
Koq = Ineq/Gfi and T ncq = h/2eV neq for current, voltage 
and time. Note that eV ncq ~ Eq/T whereas Eq. (||) is 
restricted to eV <C \/ T$cK/te,\ thus the model may lose 
its validity before or after one reaches a voltage of one 
Vneq unit, depending on which of the large parameters 
T/E c or y/T E h/E c is larger. 

The results for large te lie remarkably close to the dc 
/ — V curves of the RSJ model ^7j, except for the fi- 
nite zero-bias conductance, equal to Gn + TE,(2e/h)I neq , 
or 1 + r in the units used in the figure. Indeed, for 
small currents and voltages one may follow the analy- 
sis of Eq. (|), and find that f~E- te (dE/dip) (dcj)/dt) , 
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and G{<j>) = G N + {2e 2 NT E /h 2 T)(dE rep /d(j)) 2 ; when in- 
tegrated over a full period, our specific choice of form for 
E Tep {4>) leads to the simple 1 + r result. 



III. APPLICATIONS, LONG JUNCTIONS 

As in the RSJ model | p6| , one may use Eqs. (0) and 
(0) to study a variety of more complicated physical situa- 
tions, of which three will be considered here: the response 
of the junction to rf fields, the effects of noise, and the 
case of long junctions. 

In the presence of a high-frequency external periodic 
perturbation with frequency to, the dynamics of the junc- 
tion may phase-lock with the external drive, produc- 
ing special features called Shapiro steps at the voltage 
Vd c = huj/2e or rational multiples thereof. As in the 
previous section, the derivation for a voltage-biased sit- 
uation is simpler and may be done analytically in a 
variety of limits, but in many experiments a current- 
biased description is more appropriate We therefore 
take I = /dc + lacCosujt, and integrate Eq. (ffl) numeri- 
cally. A sample of the results is given in Fig.^(a). The 
dominance of the second harmonic, u> — 2(2e/h)V, is 
clearly displayed by the strong half-integer Shapiro step 
at V = h{fru)/2e). Note that when the energy-relaxation 
is fast, ujte <C 1 and the ac drive is weak, one can ap- 
ply Eq. (|^) and Shapiro steps cannot occur — the ratio 
of the dc components of the current and the voltage is 
in this case ^ f * G(<j>) d<j> by definition p!o| . In prac- 
tice, because of deviations from the idealized limit, the 
NEQ Shapiro steps grow linearly with u> (or Vd c ) at small 
frequencies. This was indeed observed experimentally in 
Ref. §. 

In order to model noise, one adds fluctuating terms, 
If and Jp, to the derivatives of F in Eq. (0), with 
(7 F (t)J F (0)) = 2T,pG N S(t) and similarly (J F (t)J F (0)) cx 
T 8(t). For thermal noise, = T; high-frequency exter- 
nal current noise gives > T. Again, the model may 
be integrated numerically, and Fig. ||(b) displays some 
representative results. 

The thermal case with / — > is particularly simple, 
because one may integrate out dynamical variables. If 
/ is integrated out, which is also justified for fast re- 
laxation (te — > 0), then by definition we reproduce the 
RSJ model, which in the present high-temperature case 
is trivial because I c = 0. It is more interesting to take 
te — * oo, and treat <j> as the "fast variable" , see Fig. ||(a). 
For <C E neq , the probability density of / has two peaks 
[Fig. ||(b)], and the "effective potential" F(f) has valleys 
of depth roughly ~ [Fig. ||(c)]. In the case of thermal 
noise, T = and this only leads to a complicated de- 
scription of the previous trivial case, but if T <C the 
rate of phase-slip is drastically reduced: / is attracted to 
±1, and <fi remains near the extrema of E. This explains 
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(a) An added ac current bias (with amplitudes 
I ac , frequency ui = 2eV ncq /h 1 and no noise) produces 
Shapiro steps in the dc I — V curves (here r = 25). The 
step at V = 0.5 is due to supercurrents at twice the 
Josephson frequency, (b) The effect of external noise 
(indicated by Ta) on the dc I — V curves with r = 25 
(the waviness is due to the finite numerical integration 
time). Inset: on a much expanded scale the — 0.25 
curve is seen to cross the one for Ta, = 0, and appears 



to have a "critical current" of ~ 0.1/ : 
counteracts this (T > 0; dashed lines). 



ncq- 



Thermal noise 



the very large slope displayed in the inset of Fig. [| for 
T < < -Encq — surprisingly, noise can enhance the 
zero-bias conductance here, as opposed to its effect in the 
RSJ model. Upon further expansion of the abcissa, one 
may observe that the slope never increases beyond ~ t 2 , 
despite the fact that the present discussion implies that 
this zero-bias conductance is exponentially large. This 
is due to the fact that the fluctuations in <f> induce small 
fluctuations in / (even when T = 0), which are of order 
1/te and have been ignored above. 

In junctions with a length L y much longer than the de- 
phasing length yjDrip (in the experiments of Refs. 
L y ~ 100/xm), it is reasonable to expect the dynamical 
variables / and to depend on the position y. In fact, 
a spatial dependence of <p may have a prohibitive price 
in energy, as it is associated with magnetic fields. To 
analyze this situation, the free-energy of Eq. (H) is gen- 
eralized to: 



Jo v dy dyJ 2e L y 



(8) 



where the gauge-invariant phase difference, <f>(y), and 
the local occupation probability of low-lying states, f(y), 
have become functional variables, and 



n 



, /; 0', /') = a (/ - £(</>)) 2 + /^' 2 + 7 /' 2 



(9) 
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FIG. |f (a) The ranges of noise-induced (^-fluctuations 
for fixed / (arrows), with contours at / = E (dotted 
line) and F(6, /) = 0.04E ncq = T (dashed lines; 1 = 0). 
(b) The probability density in /, for te — > oo, inferred 
from (a) is p(f) oc / d6 exp(— F/T$). (c) The effective 
potential, F(f) = — T^logp, has two dips which may 
"trap" the system, (d) The effective free energy of a long 
junction, Eq. (pi]), (c) The resulting spatial variation of 
an applied magnetic field. 



Eq. (m) is then replaced by 
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Here a = E ncq /L v , f3 = (1 /8vr) (he/ '2e) 2 (L z / L x ) is the 
energy-density of the magnetic fields associated with the 
spatial dependence of the gauge-invariant phase, H oc 
(86/ dy) (Ref. (28)), and 7 = ar^D v is introduced to 
give diffusion of / in the y direction. For simplicity, we 
set |29] the coefficient D y equal to D. 

We now turn our attention to some interesting effects 
generated by thermal noise in long junctions. Fluctua- 
tions in the occupations, i.e. the space-dependent gen- 
eralization of Af ~ TV -1 / 2 , produce fluctuations in the 
supercurrents 5F/56 (as in the case of Johnson-Nyquist 
noise, these are related by a fluctuation-dissipation the- 
orem to the zero-bias conductance, which in turn is en- 
hanced as discussed in the previous section pOjj ). The 
magnetic fields generated by these fluctuating supercur- 
rents, and the accompanying y-dependent fluctuations 
of 6, generate a y-dependence also in E. As a result, / 
relaxes to a value of E which is spatially averaged over a 
length-scale \/Dte, and can deviate from its local value, 
giving rise to a finite energy density. Indeed, after rewrit- 
ing the y integral in Eq. (p|) as a sum over wavenumbers, 
L y ^2 ql one may integrate out / one wavenumber at a 
time, which leads to a replacement of the first and last 
terms of Eq. (g) by a (l - (1 + Dr^q 2 )' 1 ) \E q \ 2 , (this is 
seen by completing the square in the exponent). The next 
step is to assume that T^JDt^ (the experimentally 
relevant regime), and to introduce a cutoff A such that 
T//3 «A« (Dte)- 1 / 2 . The fluctuations of 6 over re- 



gions < A -1 , are small because of the large value of f3 
assumed, and one may therefore linearize E(6) within 
each such region, i.e. take E q ~ (dE/d6)6 q and replace 
L y J2q by (l/27r) Jdydq. One may now integrate over 
the components of 6 with wavenumbers > A, taking care 
to keep and re-exponentiate the coefficients c in the inte- 
grals of the type J d6 q exp(—c\6 q \ 2 ), as they now contain 
an essential ^-dependence. For 7 -C j3, the resulting free 
energy is given, up to a constant, by |3l| 



F(6, 



T 



2^D^ 2/3 V , 



f dE_Y 



136" 



(11) 



where T is defined as in Eq. (g), and 6 is assumed 
smooth on the scale A -1 . The system thus has a non- 
vanishing energy density even when 6' = 0, and differs 
from a conventional long Josephson junction only in the 
replacement of — cos 6 in the usual description by a in- 
dependence with two equal-valued minima per period, 
see Fig. |(d). 

In this model, an external magnetic field (H z ) produces 
Josephson-like vortices containing half a flux quantum 
each, hc/Ae, see Fig. |](e). This may explain the mag- 
netic period-halving effects of Ref. § . The "Josephson- 
like penetration depth", defined by A 2 eq = 2y 'Dte/3 2 / '7T 
"1), is in between the exponential decay-length of 



(Ref. 

the magnetic field in regions with 6 ~ (mod 2tt) and 

the smaller one for 6 ~ n regions. The latter two may 

A 2 

be obtained by multiplying A 2 eq by (\E' )" ^ 1, where 
the primes denote differentiation with respect to 6. True 
period-halving occurs only when E ~ cos 6, e.g. when 
h/r v > Eq, and NEQ effects are weak. 

If L y -C A neq , then 6(y) ~ const. The junction prefers 
both 6 ~ and 6 ~ tt values, with intervening energy 
barriers of 2n(2e/ he) 2 ^/Dt e TE neq (L x / L z ) [This is just 



the coefficient in Eq. ( |11| ) multiplied by L y — the max- 
imal value of E' is here equal to 1]. In dc measure- 
ments which are not phase-sensitive, this behaves as a 
Josephson coupling energy, which persists to high tem- 
peratures, T S> Eq, and is roughly ~ (D/10 cL z )Ej(0) 
(Ref. p^| ). This estimate is obtained by assuming that 
\ADte is equal to a few times L x . If indeed \Jte does not 
depend strongly on temperature, one obtains a possible 
mechanism for the unexplained coupling energy observed 
in the experiments of Ref. M . 

The behavior of 6(y) in this last case is analogous to 
that of a long polymer diffusing in a disordered gel |33| : 
the polymer may be pinned by "attractive" disordered re- 
gions, leading to exponential suppression of diffusion — 
whenever one part of the polymer begins diffusing away 
from the pinning center, it is pulled back by forces due 
to the thermal motion of its other parts. If the poly- 
mer were replaced by point-like monomers, the diffu- 
sion of its constituents would not be hindered, and the 
"attractive" nature of such regions would not be notice- 
able. Likewise, if the junction of length L y <C A ncq is 



G 



replaced by much shorter junctions of length < ^/Dte, 
the high-temperature Josephson coupling discussed here 
would disappear, and only the enhanced conductance dis- 
cussed in the previous section (see Fig. m would remain. 



IV. CLEAN JUNCTIONS 

The starting point for the theory developed here was 
the fact that the spectrum of an S-N-S junction depends 
on the phase 0, but only the case of a diffusive junction 
was made explicit. Motivated by the experiments 
we will next discuss a specific clean (or qusiballistic) sys- 
tem, with N±L X 3> l e i ^ L x . Consider first the density 
of states for a rectangular geometry in the ideally clean 
case. The spectrum may be obtained for each transverse 
channel separately, by using Bohr-Sommerfeld quanti- 
zation and introducing the ^-dependence through the 
phase-shifts associated with the Andreev reflections at 
the two N-S interfaces |34]. This gives 



N± oo 



'M)=££ £ 5[e-^((2m-l),± 



±0 i—1 m— — oo 



(12) 



where Vi is the component of the Fermi velocity in the x 
direction for the ith transverse channel. We will assume 
for simplicity that the junction is three-dimensional, i.e. 
that there are two substantial transverse directions. The 



sum over transverse channels, ^ 



N± 



may then be re- 



p, where 



placed by the integration N± J Q " F . . . 2v± dv± / v 
v i = \J v f ~ v \ ■ This has the effect of replacing every 



(5-function in Eq. (12) by a "sawtooth": 



^££#e( £m , ± 



(13) 



c m,± 



where e mi ± = \Eq ({2m — l)n±(f>), the Thouless energy 
is Eq — hvp/L x , and is the Heavyside step-function. 
This density of states rises linearly at small energies , 
see Fig. ||, as opposed to the mini-gap of the diffusive 
case [for a two-dimensional junction, each term in the 
sum would be divided by 2yl — e/e m ,±, introducing di- 
vergencies in v(t,4>)\ as these divergencies are integrable, 
the NEQ results would still be qualitatively similar]. 

Weak elastic scattering will result in avoided crossings 
in the spectrum. The disorder-averaged excitation en- 
ergy of the nth state may then be introduced as above, 
n = J Q n deu(e, (f>). These are displayed in Fig. |, which 
is to be compared with the diffusive case of Fig. n , We are 
assuming here that Tl/t c \ <C Eq and therefore the elas- 
tic scattering does not appreciably smear the features 
in i/(e, </>). Note that the parameter r e \ is nominal and 
describes the scattering rate in the absence of proxim- 
ity effects; the proximity effects change the scattering 
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FIG. g: (a) Energies E n {4>) of representative Andreev 
bound states, for a (three-dimensional) quasiballistic S- 
N-S junction, as in Fig. [|. Here "quasiballistic" means 
that there is enough disorder to produce avoided cross- 
ings (not resolved), but not enough to modify the den- 
sity of states appreciably from that of the super-clean 
limit, (b) The density of states of the junction, sam- 
pled at <j> = 0.37T. The NEQ supercucrrents in the limits 
of small and large {2e/Ti)T^y are plotted in (c) and (d) 
respectively. The model using a simple representative 
curve (dashed lines) is compared to the result of the full 
sum over the different E n (cj>) levels (full lines). 



rates, which are given by the product of v(e, (f) with the 
variance of the elastic-scattering matrix-elements. Thus 
the typical change in e suffered by an individual level 
due to the disorder may be substantially increased or de- 
creased relative to Ti/t c \. This is particularly relevant in 
the region near <p = ir and e = 0, where all of the ./Vl 
states associated with the lowest-lying term in Eq. (|TJ) 
become crammed into a very small range of energy. We 
will ignore such complications here — their resolution 
(e.g. the size of the region in e and <fr which is affected) 
depends on the value of Eqt c i/H which we are taking to 
be infinite. In fact, one should also note that apart from 
assuming A ^> e, the Andreev approximation used in de- 
riving Eq. (|l2|) is only valid when the kinetic enery in 
the x direction, ^m*vf, is much larger than the gap A 
(here in* is the effective mass of the electrons). In prac- 
tice, a fraction A/_Ep of the states do not have a large 
enough Vi\ although they are not strongly affected by 
Andreev reflections ("glancing angle") and show much 
weaker proximity effects |3~i|| , they are ignored here by 
taking Ep/A to infinity. 
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We next consider dynamic NEQ effects in this system. 
For brevity only the results in the simple limits of small 
and large voltages relative to the inelastic scatterring 
rate, Eqs. (||) and (Q), are plotted in Fig. ||(c) and (d). 
Note that at moderatley large voltages, (2e/K)V ~ l/r e i, 
the avoided crossings may "no longer be avoided" , and 
the ballistic structure of the spectrum may substantially 
alter the analysis. In other words, the occupations might, 
under sufficiently strong drive, become dependent not 
only on energy but also on the angle of propagation rel- 
ative to the x direction (i.e. the transverse channel index 
i). As explained in Sec. II, the use of the relaxation- 
time approximation, Eq. (|^), is also expected to fail at 
large voltages, when (2e/K)V ~ y/EcJhrE. Either one 
of these restrictions may be more stringent for a partic- 
ular application, because te is usually much larger than 
r e i, which itslef is 3> Ti/Eq in the clean case. 

As in the case of the diffusive junction, it is pro- 
posed here that the complicated energy-dependence of 
the spectrum can be ignored to a first approximation, 
with the phase-dependence and the NEQ occupation fol- 
lowed only for a single representative level. The individ- 
ual curves E n ((f>) in Fig. ^|have discontinuous derivatives 
(because v is discontinuous), and are not represented well 
by any single one of them. Instead we choose the form 

E((f>) = 1 — 2^/| sin(i0)| — i| sin(0)|, which reproduces 

the qualitative behavior of G(</>) near both = (where 
it vanishes linearly) and <f> — 7r (where it approaches a 
constant). This form has no adjustable parameters, and 
gives errors of order 30%, as shown in Fig. || [the value of 
E ncq used was 2.Q{Eq/8T)\. These errors could of course 
be reduced by using a more complicated form of E(4>), 
but it turns out that the contribution of the levels near 
e = ttEq end <\> — is difficult to model with a single 
representative curve (to correct this, one could perhaps 
introduce a ^-dependent number of represented levels N, 
which would vanish linearly at small 0). Despite these 
differences in E(4>), and even the possible need to follow 
the occupations of many levels instead of using a sin- 
gle representative, we expect that the results plotted in 
Figs. |^ and ^ would not be dramatically altered if the 
calculation were redone for a clean junction. The results 
for long junctions, Fig. would obviously be more sen- 
sitive to the qualitative behavior near = and — n. 
However, for any realistic system there will be additional 
effects which would "soften" this behavior (see below). 

V. DISCUSSION 

In summary, it was shown here that in wide S-N-S 
junctions, significant deviations of the electronic occu- 
pations from equilibrium may be caused not only by an 
external dc or ac drive, but also by external noise, or 
even by space-dependent thermal noise in the case of 
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FIG. |: Sketch of the different versions of E((f>) suggested 
for a diffusive junction (full line, Sec. II), a clean junc- 
tion (dashed line, Sec IV), and the simple cosine used in 
modelling the experiment (dotted line, Ref. 0), shifted 
for clarity. The "bottom of the valley" shown in the inset 
of Fig. U and in Fig. ||(a) should be modified accordingly. 

long junctions. Such non-equilibrium occupations have 
structure in the excitation energy e on the scale of the 
Thouless energy Eq , and therefore have "ensemble coher- 
ence lengths" large enough to span the width of the junc- 
tion L x . For this reason, even when the temperature T is 
high, and the equilibrium "ensemble coherence length" , 
£nj is small enough to give an exponential supression of 
the equilibrium Josephson effect, supercurrents may flow 
in the junction. These supercurrents lead to Josephson 
like characteristics of the junction, and are charactrized 
by an approximate period-halving and a weak tempera- 
ture dependence. 

The model developed here to describe these effects, 
defined by Eq. (^) and Eq. (0V has been compared with 
experimental results in Ref. ||], for both dc and ac ef- 
fects, [Fig. H and Fig. ||(a) respectively]. It was found 
there that the high harmonics are very weak, and the 
model with E = cos (see Fig. 0) describres both types 
of effects well with the same choice of parameters, I noq 
and te (the latter was further corroborated by estimating 
it from the sub-gap structure). This may serve as addi- 
tional motivation for representing all of the states in the 
spectrum by a single curve -E(0), as was done here — the 
fine details which are being left out by this procedure are 
difficult to observe experimentally. The most important 
factor which contributes to the "softening" of the phase- 
dependence, compared to the clean system (l c \ ^> L x ) 
considered in the previous section, is probably the elas- 
tic scattering, l c \ ~ L Xl in the experiment. 

Theoretically, there are many refinements of the model 
which remain to be pursued, of which three are listed 
here: (a) No distinction has been made here between 
electron-electron and electron-phonon inelastic scatter- 
ing; instead all mechanisms which smooth the occupation 
functions on the scale of Eq have been lumped into one 
parameter, te- Although this may be adequate, the rela- 
tionship of this parameter to r c _ c and r c _ p h needs to be 
examined more carefully, especially as it will determine 
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its temperature dependence, (b) The low-voltage limita- 
tion of the present model, eV <C yj Eq W/te, may perhaps 
be relaxed to eV <C Eq without significantly complicat- 
ing the model, by adding a phenomelogical term involv- 
ing £>e oc V 2 to the evolution equation of /, Eq. (||). 
This is important because experimental I—V curves of- 
ten exhibit interesting features in this range of voltages 
M. (c) The lateral diffusion of / in long junctions has 
been treated here in a phenomelogical manner, by setting 
D y = D; a more sophisticated treatment of this may be 
in order (cf. Ref. 

Clearly there is a richness of phenomena in such junc- 
tions which remains to be explored experimentally as 
well. The most dramatic challenge, perhaps, is to observe 
the predicted "half flux quanta" of Fig. ||(e) directly, e.g. 
by pinning them and using magnetic microscopy tech- 
niques similar to those used in observing half flux quanta 



in junctions made of high-T c materials 35 
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APPENDIX: COMPARISON WITH THE 
GREEN'S FUNCTIONS APPROACH 

In order to compare our results with those obtained 
with the Green's function approach in the literature, we 
first quote the Usadel equations, in the form used by 
Zhou et al. jjj . In the normal-metal part of the diffusive 
junction of Fig. |], the equations read: 



D 
2 l 



+ ie sm( 



D 



( X ) sin 29 = ; 

(x'sm 2 ey = o, 



(14) 



in terms of the variables 9 and x, which are both complex 
and depend on position x and excitation energy e (the 
primes denote differentiation with respect to x). The 
excitation energy e here is taken to be positive (at e = 
the variables \ and 9 are real and can be thought of as 
angles). Here Ti has been set to 1, and D is the diffusion 
constant as before. The boundary conditions are 9 = ^ir 
and x = ±5^ a t the N-S interfaces, x = ±^L X . For a 
recent derivation and numerical implemetation of these 
equations, see Ref. pl[ . 

Once the Usadel equations have been solved, one has 
the Green's functions (in angular variables), and may ob- 
tain many of the properties of the junction. Specifically, 
the local density of states is given by 5fe cos 9 times its 
value in the absence of proximity effects, and the Joseph- 
son current density, j(e,(j>), is given by the imaginary 
part of x' sm2 9- The fact that this must be equivalent 
to the mesoscopic-systems approach used in Eq. (Q), i.e. 
that j(e, <t>) — {2e/K)v{dE n /d(j)) is clear if one recalls that 



the phase cj> can be represented through a vector poten- 
tial, and that a derivative of the energy operator with 
respect to the vector potential gives the current opera- 
tor. It is a non-trivial excercise to derive this relationship 
analytically from the disorder-averaged Eq. (|l4|), but it 
is easy to check it in practice once the solutions have 
been obtained, and this has in fact been applied to eval- 
uate the accuracy of the numerical scheme used to plot 
Fig. [l| Note also that the normal current 7n appears in 
the Green's function approach, e.g. in Ref. |TJ], in exactly 
the same way as in Eq. (|l|). 

We next turn to the evolution of the occupations f n . 
In order to represent both electrons and holes without 
extending e to negative energies, the Green's-functions 
approach uses two occupation functions: the symmetric 
combination fi(e,x) = 1 — f(—e,x) — f(e,x), and the 
antisymmetric combination f2(e,x) = f(—e,x) — f(e,x) 
[note that the sum over n in Eq. (|l|) may be restricted to 
positive energies by replacing /„ by —/a, because EL« = 
— E n ]. In equilibrium, one has f% = fo = tanh(e/2T) 
and /i = —etpdf2/de, where a small scalar potential 
y>(x) has been introduced. In the present work, we have 
essentially introduced the electric field through a time- 
dependent vector potential, so that f\ — throughout 
(cf. Ref. p3]). The evolution of f 2 , which is space- 
independent for eV <C Eq, is controlled by 



dh + dh dEn d± = _ 1_ 
dt de d(j> dt te 



fo 



(15) 



This is just Eq. (^|) rewritten using the variables e and t 
rather than n and t as independent (we have = 1 — 
2/„ when e = E n ). By noting that (dE n /d<j))(d<j)/dt) = 
(j/p)V, one finds that Eq. ( |l5| ) is identical to Eq. (12) of 
Ref. flifl , where a scalar potential with cp = (x/L x )V has 
been used instead of a vector potential, and therefore the 
difference between the values of /i evaluated at the two 
N-S interfaces is equal to eUd/2/de (Ref |p6| ). 

We have thus shown that Eqs. ([!]) and (|2j), which form 
the basis for the present work, are equivalent to those 
underlying more conventional treatments of NEQ effects. 
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